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Abstract. We study the production of charm mesons and other charm baryons 
from quark matter at CERN SPS and RHIC energies. Using quark coalescence 
models as hadronization mechanism, we predict particle ratios, absolute yields 
and transverse momentum spectra. 
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1. Introduction 



During the last decade J/ip production was investigated with great interest at the 
CERN Super Proton Synchrotron (SPS), and an anomalous suppression has been 
detected in Pb+Pb collisions at 158 GeV/nucleon bombarding energy jjj. This result 
is is assumed to be the strongest evidence of QGP formation (for an overview see 
Ref. Q). However, the unambiguous confirmation of the QGP formation requires 
more charm related data, e.g. absolute yields, transverse momentum spectra and the 
characteristics of charm hadron production from deconfmed matter. 

In this paper we summarize our results obtained by assuming quark coalescence 
as microscopical hadronization mechanism of the possibly produced deconfined phase 
at CERN SPS energy. The ALCOR |, |, § and MICOR |, § hadronization models 
were extended by including charm quark and the consecutive charm hadron formation. 



2. The number of charm quarks in heavy ion collision at SPS energy 

Perturbative QCD calculations predict small number of newly produced charm quark- 
antiquark pairs in the Pb+Pb collision at CERN SPS energy, N C c « 0.1. On the other 
hand statistical models obtain a much larger number, N C c ~ 8.5 If we assume that 
deconfinement appears after the formation of a hot dense string matter where strings 
become ropes in the presence of large density (characterized by an increased string 
constant, k), then quark-antiquark pair creation will be dominated by the breaking 
of these in- medium strings/ropes. The enhancement of the heavy quark-antiquark 
production, jx, can be calculated by means of the Schwinger mechanism |^|: 

P(XX) ( Tr(m x -m 2 g ) \ 

7 ^W =CXP ( * J 

If we increase the usual value of the string constant k = 1.0 GeV/fm with 10 %, 
then an enhancement factor 24 appears for charm quarks with mass m c — 1.5 GeV 
(m q — 223 MeV). 20 % increase in k would generate an enhancement factor 350 !!! 
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Thus relatively small change in the string breaking could modify dramatically the 
total charm production in heavy ion collisions. If relatively large number of charm 
quarks can be produced in the initial phase of the heavy ion collision via in-medium 
breaking of the primary strings, then after deconfinement the quark coalescence can 
become a reliable description of charm hadron formation from the plasma state. 

3. Quark coalescence 

In our description quark coalescence is the basic hadronization mechanism. It is 
assumed that the number of directly produced hadrons is given by the product of the 
number of quarks (or anti-quarks) from which those hadrons are produced, multiplied 
by coalescence coefficients Ch and by non- linear normalization coefficients b q , that 
take into account conservation of quark numbers during quark coalescence [|[ |] . The 
number of various quarks is denoted by the symbol usual for that type of particles, e.q. 
u, d, s and c denote the number of light, strange and charmed quarks, respectively. 
The following examples display the appropriate coalescence relations (see Ref. ||): 
p + (uud) = C p ■ (b u uf ■ (b d d) K+ (us) = C K ■ (b u u) ■ {hs) 

n (sss) = Cn ■ (b s sf 4> (ss) = Cp ■ (b s s) ■ (b^s) 

n c (ssc) = C^-(b s s) 2 -(b c c) D+(cd) =C D -{b c c)-{b 1 d) 

n cc (scc) =C^-(6 sS )-(6 c c) 2 D+(cs) =C s D -(b c c)-(bjs) 

n ccc (ccc) = Cg c ■ (b c c) 3 J ft) (cc) = Cj/^ ■ {b c c) ■ (he) 

The non- linear coalescence factors b u , bd, b s , b c and the bu, b^, bj, be are determined 
unambiguously from the requirement that the number of constituent quarks and anti- 
quarks do not change during the hadronization, and that all initially available quarks 
and anti-quarks have to end up in the directly produced hadrons. This constraint is 
a basic assumption in all models of quark coalescence. The correct quark counting 
yields to 8 coupled non-linear equations, expressing the conservation of the number of 
quarks and antiquarks. Without solving these non-linear equations, one can recognize 
interesting relations among the different hadron ratios. For example certain meson 
ratios play the role of step functions among different antibaryon to baryon ratios. 
This is illustrated in Figure 1 from Ref. Q , where the u and d quarks are considered 
together as q light quark. 




s qs 

Figure 1. Mesonic step factors S qs , S sc , S cq and antibaryon to baryon ratios. 
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Here the step factors are defined as: 
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These properties lead to the identity S qs ■ S sc ■ S cq = 
as a measurable relation between these mesonic ratios: 
These relations are valid at both CERN SPS and RHIC energies, if quark coalescence 
describes the hadronization of the produced deconfined phase. 



1, which can be rewritten 
(D s /Ds)/(D/D) = K/K. 



4. Predictions for charm hadron yields at SPS and RHIC 

Total hadron yields can be determined through the numerical solution of the 
appropriate 8 non-linear equations of quark number conservation. The input values 
are the newly produced quark-antiquark numbers, N u ^, N^, N S g and N C c- At CERN 
SPS energy the earlier ALCOR calculations could successfully reproduce most of the 
experimental data within errorbar (even the fl/fl ratio) Q. 

Including the charm quark into this model, one experimental data about charm 
hadron production is needed to determine the full charm production yield. We choose 
the known ratio J/ip/h~ = 1.0 ■ ||. The ALCOR model can reproduce this data 
with N C c = 3.4 ]lQ] . The obtained charm hadron yields (see Table 1) display a factor 
of 12 enhancement relative to the perturbative QCD predictions, but they reach only 
40 % of the predictions from the statistical model Q . 

We repeat our calculation at RHIC energies (see Table 1), assuming that the 
relative strange and charm quark-antiquark pair production remains the same, only the 
absolute yields are scaled with increasing energy, following the increase in the number 
of h~ (h~ was obtained from the HIJING model). Table 1 shows that detectable 
numbers of charm hadrons (D, A c , A c ) ought to be produced at RHIC. 



5. Transverse momentum spectra of the D and J/tp mesons 

The transverse momentum spectra of the produced hadrons can be calculated by 
the MICOR model §, §. Applying this model one can determine the hadronization 
temperature Thad and the transverse flow of the hadronizing quarks. At CERN 
SPS energy T HAD = 175 MeV and v T = 0.46 ± 0.05 fit the data in the Pb+Pb 
collision |3| . Including charm quarks into this model we can determine the transverse 
momentum spectra of the D and the J/ip mesons. Figure 2 displays the obtained 
values T® ff = 259 ± 19 MeV and Tjjf = 315 ± 35 MeV for the slope parameter of 
the transverse momentum spectra |j. We urge future experiments to measure these 
quantities at CERN SPS. 
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Table 1. Charm hadron production at CERN SPS and RHIC energies from the 
ALCOR model 
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Figure 2. Experimental hadronic slopes of the transverse momentum spectra in the 
PbPb collision at 158 AGeV energy from WA97 (squares), NA50 (dots) and NA44 
Collaborations (triangulars), see Ref. §. Open squares with errorbar indicate the 
MICOR results on D and Jji\> mesons. 
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